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ABSTRACT 

Benchmark aerodynamic data are presented for compressible flow through a repre- 
sentative S-duct configuration. A numerical prediction of the S-duct flow field, obtained 
from a subsonic parabolized Navier-Stokes algorithm, is also shown. The experimental and 
numerical results are compared. 

Measurements of the three-dimensional velocity field, total pressures and static pressures 
were obtained at five cross-sectional planes. Aerodynamic data were gathered with calibrated 
pneumatic probes. Surface static pressure and surface flow visualization data were also 
acquired. All reported tests were conducted with an inlet centerline Mach number of 0.6. 
The Reynolds number, based on the inlet centerline velocity and duct inlet diameter, was 
2.6 x 10 6 . Thin inlet turbulent boundary layers existed. 

The collected data should be beneficial to aircraft inlet designers and the measurements 
are suitable for the validation of computational codes. The results show that a region of 
stream wise flow separation occurred within the duct. Details about the separated flow region, 
including mechanisms which drive this complicated flow phenomenon, are discussed. Results 
also indicate that the duct curvature induces strong pressure driven secondary flows. The 
cross flows evolve into counter-rotating vortices. These vortices convect low momentum 
fluid of the boundary layer toward the center of the duct, degrading both the uniformity and 
magnitude of the total pressure profile. 
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CHAPTER I 
INTRODUCTION 

Aircraft propulsion systems often use S-ducts. A primary purpose of an S-duct is to 
convey air flow from the wing or fuselage intake to the engine compressor. Examples of 
commercial aircraft with S-ducts include the Boeing 727 and the Lockheed Tristar L- 1011. 
Amongst military aircraft, both the General Dynamics F-16 and the McDonnell-Douglas 
F-18 use S-shaped ducts. 

Diffusing S-ducts have centerline curvature and cross-sectional area increase. Curvature 
of the centerline or changes in the duct cross-sectional shape give rise to streamline 
curvature. Cross-stream pressure gradients, resulting from streamline curvature, can produce 
significant secondary flows. Additionally, the adverse streamwise pressure gradient, caused 
by increasing cross-sectional area, can lead to flow separation. 

Well designed diffusing S-ducts should efficiently decelerate the compressible incoming 
flow. Moreover, to achieve appropriate performance, the S-duct must also incur minimal total 
pressure losses and deliver nearly uniform flow with small transverse velocity components 
at the engine compressor entrance. Often the aircraft designer faces a difficult dilemma. 
Size and weight restrictions encourage the use of shorter S-ducts. However, this results in 
greater streamline curvature, larger adverse pressure gradients and the risk of unacceptable 
duct performance. 

The ability to quickly predict complicated flow phenomena without experimental testing 
is obviously desirable for designers. The highly three-dimensional flow in a diffusing S- 
duct presents a substantial challenge to computational fluid dynamic (CFD) flow solvers. 
Unfortunately, current numerical algorithms normally require a large amount of computer 
time to solve the Navier-Stokes equations, which are required for complex flows. Moreover, 
only a limited amount of experimental benchmark data exists for the validation of CFD 
codes applied to compressible S-duct flows. 

This report documents the important details of a concurrent experimental and computa- 
tional study of the compressible flow through a diffusing S-duct. Experimentally, internal 
flow measurements of the three-dimensional velocity field, total pressures and static pressures 
were obtained in five cross-stream measurement planes. Surface static pressures and surface 


2 


flow visualization were also used to reveal flow field characteristics. Numerically, the flow 
field was modeled using an efficient subsonic parabolized Navier-Stokes (PNS) algorithm. 

The emphasis of this report is placed on the experimental measurements. The objective 
of the experimental study is to investigate the development of the flow field and to provide 
a comprehensive benchmark data set for the compressible flow through a representative 
diffusing S-duct. Numerical calculations are presented concurrently with the measured 
results to help explain the development of the flow. The computational results also reaffirm 
that PNS models quickly and efficiently produce qualitatively correct results of complex 
internal flow fields. 

This study could be beneficial to inlet designers. The numerical results show that flow 
field trends can be predicted using the simplified parabolized model. The measured data 
are suitable for CFD code validation. Data for S -ducts containing active or passive flow 
management devices intended to improve duct performance could be compared to these 
baseline data. 
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CHAPTER II 

A REVIEW OF PRIOR INVESTIGATIONS 

Previous experimental and computational studies have been conducted to determine the 
character of the fluid motion within S-ducts. From the available literature, selected studies 
similar to the one presented in this thesis are reviewed. Attention will be drawn to both 
experimental investigations and computational studies. 

The design and implementation of earlier aircraft propulsion S-duct components relied 
heavily on test-stand data, because, in the past, these flow fields were too complex for 
numerical modeling. The development of the Boeing 727 center engine inlet duct was 
documented by both Bauermeister et al. [1] and the Boeing Commercial Aircraft Company 
[2]. Preliminary duct designs of the 727 inlet contained regions of reversed flow and 
had unacceptable exit pressure distributions. Center engine surging was also experienced 
during initial flight testing. A detailed experimental study revealed how installation of 
co-rotating vortex generators within the duct improved exit pressure uniformity, prevented 
flow separation and provided surgefree operation. In a separate study, the performance of 
several S-duct scoop diffusers used in conjunction with high-speed turboprop-fan engines 
was examined by Little and Trimboli [3]. These S-ducts had severe centerline curvature 
and underwent radical changes in cross-sectional shape. The study suggested the existence 
of a complicated three dimensional duct flow and provided valuable performance data for a 
range of applications. Both of these studies are prime examples of the extensive amount of 
experimental research needed to design efficient S-ducts used for engine inlets. 

To aid in inlet design, experimental studies have also been performed on S-ducts 
which represent simplified models of engine inlets [4-9]. These studies were conducted 
to qualitatively examine S-duct flow fields and provided better understanding of the flow 
physics. These studies found counter-rotating vortices within S-ducts. The vortices convected 
low momentum fluid toward the center of the duct. A highly nonuniform exit velocity field 
existed. Rowe [4] theoretically proved that the development of the vortices was an inviscid 
phenomenon, provided that a nonuniform inlet velocity profile (typical of a boundary layer 
profile) existed. Taylor et al. [10,11,8] showed that counter-rotating vortices developed 
in ducts having circular or square cross sections. In addition, these vortices were present 
regardless of whether the inlet boundary layer was turbulent or laminar. S-ducts having large 
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incident inlet flow angles (Guo et al. [7,12]) were absent of counter-rotating vortices and 
instead the entire exit flow rotated, similar to solid body rotation, about the centerline. 

Unfortunately, the majority of the studies involving simplified duct models were 
conducted with incompressible flow and/or thick boundary layers. Furthermore, only a 
few involved diffusing ducts. Most of the studies did not have streamwise separated flow, 
because most of the ducts were non-diffusing and had only mild centerline curvature. 
Streamwise separated flow did occur in the ducts studied by Bansod and Bradshaw [5]. 
However, the region of separation was small and had no significant effect on the flow field. 
A large separated flow region was present in a study by Sullivan et al. [13], but only 
surface flow visualization and surface static pressure data were obtained. Further studies 
involving representative diffusing S-ducts, having streamwise separation, compressible inlet 
flow conditions and thin inlet turbulent boundary layers, were needed to help interpret the 
flow physics associated with S-ducts used for aeronautical applications. 

Studies on simplified duct models have also been conducted to obtain experimental 
benchmark data, which are needed to determine the effectiveness of a numerical model’s 
ability to represent the flow field. A considerable amount of benchmark data exists 
for incompressible flow through S-ducts [7,10,11,8]. The only detailed measurements of 
compressible flow through a diffusing S-duct were acquired by Vakili et al. [9]. In this 
study, surface flow visualization suggested the existence of three-dimensional separated 
flow. Aerodynamic data were obtained with a non-calibrated five-port cone probe. These 
measurements showed that the exit flow was nonuniform. 

The importance of computational fluid dynamics (CFD) as a design tool continues to 
be established. CFD internal flow research includes fully elliptical Navier-Stokes (FNS) and 
parabolized Navier-Stokes (PNS) models. Past efforts were focused on PNS methods because 
of their efficiency and small storage requirements. Computations involving many different 
duct geometries have been accomplished by a NASA Lewis Research Team using a PNS 
algorithm designated PEPSIG [14—17]. Results showed, in general, acceptable agreement 
between numerical and experimental data, when experimental data were available. However, 
the computations underpredicted the extent of exit flow nonuniformity. The location and 
extent of separation regions were also incorrectly predicted. A review of the literature 
suggests that the turbulence model was a major cause of this discrepancy. 
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With the availability of faster and larger memory computers, solving the complete Navier- 
Stokes equations has become possible. Unfortunately, only a few calculations involving 
S-ducts have been reported. Smith et al. [18] and Harloff et al. [19,20] have prepared 
reports documenting the flows in non-diffusing and diffusing circular cross-section S-ducts. 
A finite-volume FNS scheme and Baldwin-Lomax turbulence model were used to predict 
the flow field. Effects of different mesh structures on the flow field were examined. The 
reports indicate acceptable qualitative agreement with Vakili’s et al. [9] experimental data. 
Nevertheless, within the diffusing S-duct, the accurate location and length of separation 
and the extent of exit flow nonuniformity were not correctly predicted. The discrepancies 
between computational and experimental data were possibly due to the turbulence model, 
the mesh density or undocumented disturbances in the experimental measurements. These 
studies suggested that a k-t turbulence model should be incorporated into the FNS algorithm, 
a finer mesh should be used and additional comprehensive experimental data should be 
obtained for compressible diffusing S-duct flows. 

The material in this thesis is a partial response to the suggestions of Smith etal. [18] and 
Harloff etal. [19,20]. Benchmark data have been obtained for the compressible flow through 
a diffusing S-duct. These data can be used to validate the computational results obtained 
from the FNS algorithm which currently incorporates an algebraic or k-t turbulence model. 
In addition, PNS results reported in this thesis enable interested individuals to compare PNS 
and FNS prediction capabilities, advantages and disadvantages. Finally, the combination of 
the detailed experimental and computational data presented together enhances a designers 
ability to interpret the complicated flow physics within a diffusing S-duct. 
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CHAPTER m 

EXPERIMENTAL FACILITIES AND PROCEDURES 

One of the goals of this study was to provide a detailed set of test flow measurements. 
The calibration facility, test facility, S-duct geometry, instrumentation and measurement 
techniques used in the experimental program are described in this chapter. 

Probe Calibration Facility 

A new facility was constructed for the calibration of aerodynamic probes in the Internal 
Fluids Mechanics Facility at the NASA Lewis Research Center. The facility was designed for 
convenient access and with it the accurate calibration of measurement probes was possible. 
The facility consisted of a flow conditioning section, calibration point and an exhaust diverter. 
A schematic of the facility is shown in Figure III.l. 



Figure IHl Schematic of the probe calibration facility 
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Flow conditioning section 

The flow conditioning section contained a series of filters, meshes and screens which 
provided a uniform ffeestream flow field. Pressurized shop air was used as the flow driving 
mechanism. Mach numbers up to 0.85 were attainable. The flow rate was controlled by a gate 
value. A Cuno filter retained foreign particles from continuing downstream. A perforated 
steel cone mixed the flow. A two inch wide honeycomb mesh reduced large eddies into small 
eddies. A fine mesh screen eliminated any local nonuniformities created by the honeycomb 
mesh. Plenum total pressure was measured by a Pitot probe positioned downstream of the 
screen. Therefore, only minimal losses occurred between the measured total pressure point 
and the calibration point. An exit contraction nozzle uniformly accelerated the flow to the 
calibration point. The diameter of the contraction nozzle at the exit was nominally 5 cm. 

Calibration point 

The calibration point was located 5 cm away from the contraction exit and in the center 
of the exhaust flow. The facility allowed rotations in the pitch and yaw directions to occur 
without repositioning the calibration point. Pitching motion was electronically controlled 
by an L.C. Smith actuator. Yawing motion was manually adjusted by the operator using 
a Klinger micro-angle-control (accurate to within 0° 01’). A photograph of pitch and yaw 
mechanisms assembled on the calibration rig is shown in Figure III.2. 

Experimental Test Facility 

All experiments were performed at the NASA Lewis Research Center in the Internal 
Fluid Mechanics Facility. This facility was designed to support tests of a variety of internal 
flow configurations. Complete details of the facility are described by Porro et al. [21] For 
this experiment, atmospheric test-cell air was drawn through a large settling chamber, passed 
through the test section and dumped into a large exhaust plenum. The use of test cell air 
inhibited independent variation of Mach numbers and Reynolds numbers. A schematic of 
the facility is shown in Figure 113.3. 

Settling chamber 

The settling chamber conditioned the incoming flow in the following way. Air was 
drawn into the chamber through a large bell-mouth opening. A perforated spreader cone 


Figure III.2 Pitch and yaw mechanisms assembled on the calibration rig 


Settling Chamber 


Test Section 


Exhaust Section 


* 





Figure III. 3 Schematic of the Internal Fluid Mechanics Facility 
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mixed the inlet flow. A coarse mesh conditioning screen reduced mean flow nonuniformities. 
A honeycomb- screen combination removed large scale turbulence fluctuations. A seamless 
contraction section uniformly accelerated the flow from the settling chamber. An area 
contraction ratio of 59 to 1 ensured a low turbulence intensity flow and nearly uniform flow 
at the test section entrance. 

Test section 

The test section for this experiment consisted of the diffusing S-duct and two constant area 
duct extensions. The first extension (10.21 cm diameter) served as the interface between 
the contraction exit and the S-duct entrance. The second extension (12.57 cm diameter) 
conveyed the flow from the S-duct to the exhaust region. The second extension was able 
to rotate about its centerline, when the facility was shut down. Each extension was 76.2 
cm long and had the same internal surface finish as the S-duct. S-duct details are review in 
sections S-Duct Geometry and Instrumentation and Measurement Techniques. 

Exhaust section 

The exhaust section contained a circular cross-section pipe, a mass flow plug and a 
sub-atmospheric plenum. The pipe was 244 cm long and had a diameter of 12.70 cm. 
The finish was unpolished. The pipe housed the adjustable mass flow plug and also 
assured no downstream influence of the exhaust plenum on the test section. The mass 
flow plug controlled the airflow through the entire facility. Mach numbers up to 0.95 within 
the test section were possible with this arrangement. Flow was simply dumped into the 
sub- atmospheric pressure exhaust plenum, which had a 121.92 cm diameter. 


S-Duct Geometry 

The diffusing S-duct was intentionally designed to incorporate as many of the complex 
three dimensional flows, including the possibility of unsteady streamwise separation, 
associated with similar configurations. The duct was designed at the NASA Lewis Research 
Center. A half shell representation of the diffusing S-duct is shown in Figure I II .4. 
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Figure III.4 Half shell representation of the diffusing S-duct 

The duct centerline is defined by two circular arcs with identical radii, R, of 102.1 cm and 
subtended angles, 0 max /2, of 30°. Both arcs lie in the x-z plane. The centerline coordinates, 
given by Equations III. 1 , III.2 and ID.3, are indicated by the dashed line in Figure III.4. All 
cross sections perpendicular to the centerline are circular. The duct inlet radius, r\, is 10.21 
cm. The duct exit radius, r 2 , is 12.57 cm. This produces an area ratio, Ao/Ai, of 1.52. The 
variation of the duct radius as a function of the angle 9 is given by Equation III.4. The duct 
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is larger than, but geometrically similar to, the duct tested by Vakili et al. [9]. 
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When discussing results, axial position refers to the distance to cross-stream planes 
normalized by the inlet diameter, measured along the duct centerline, beginning at the start 
of curvature. Positions within any cross-stream plane are specified by the polar angle <j), 
measured from the positive z-axis, and the radial distance r from the centerline. Finally, the 
S-duct is fixed to the coordinate system shown in Figure IH4. 

The duct was milled from two separate blocks of aluminium and had a final tolerance 
of ±0.0127 cm. After milling, the two halves were mated together and sanded using 120 
grit Emory cloth until all machining imperfections were removed. The surface was then 
polished. The interior split line, located on the vertical centerline plane, was flush to touch 
and invisible to sight. 


Instrumentation and Measurement Techniques 

The flow field was investigated with several different techniques. Flow visualization 
and wall static pressures were obtained on the S-duct surface. Three- and five-hole pressure 
probes were used to determine total pressures, static pressures and the velocity field. Each 
measurement technique is briefly explained within this section. 

Surface flow visualization 

Surface flow fields can be visualized by applying fluorescent oil to a duct wall and then 
observing streakline patterns which develop. The method is especially advantages to use 
when nearly steady state flow conditions exist. In general, regions of streamwise reversed 
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flow and boundary layer cross flows can be identified with this technique. The location of 
flow separation and appreciable cross-stream pressure gradients can then be deduced. 

The assumption made with surface oil flow visualization is that the oil streaklines follow 
the flow streamlines near the surface of the duct. Squire [22] concluded that oil streaklines 
follow the boundary layer surface streamlines except near separation. During separation, 
either in compressible or incompressible flows, the oil tends to form an envelope slightly 
upstream of the true separation region. However, this premature indication of flow separation 
is less marked for turbulent than for laminar boundary layers. Squire [22] also determined 
that the effect of the oil on the boundary layer flow is negligible in most practical cases. 
Reichert et al. [23] have had success using a variation of this technique to visualize subsonic 
flow in a circular-to-rectangular transition duct. Davis et al. [24] have also successfully 
used a thin-oil-film method to investigate near wall flow behavior in hypermixing nozzles. 

For this study, a mixture of fluorescent commercial powder dye (Rocket Red : Day-Glow 
Color Corporation) and petroleum-base lubricating oil (STP 140 wt.) was prepared. Slight 
thinning of this mixture with mineral oil was needed depending on the amount of powder 
dye used. The mixture was applied to the surface in small dots. This small-dot application 
method (proposed by Wendt [25]) appeared to be an improvement to the film technique 
utilized by Reichert et al. [23]. After oil application, the duct was quickly assembled for 
testing. Surface oil streakline patterns developed when air was drawn through the duct. 
The patterns corresponded to surface flow direction and shear intensity. The duration of 
the test, at the desired flow rate, was ten minutes. Following rapid shut down, the duct 
was disassembled in order to observe the streaklines under the illumination of ultraviolet 
light. The patterns were photographed with a 35mm camera using 400 ASA color film. An 
UV filter was used to reduce the glare produced by the reflection of the ultraviolet light 
from the surface. 

Pressure measurements 

The test section was instrumented to obtain mean pressure measurements. Total and 
static pressures were measured in cross-sectional planes with three- and five-hole pneumatic 
probes. Surface static pressures were acquired with wall taps. The pressure instrumentation 
and measurement system are described in this section. 
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Table in.l Experimental measurement plane information 


Plane 

A 

B 

C 

D 

E 

Location, s/d 

-0.50 

0.96 

2.97 

4.01 

5.73 

Radial Traverses 

10 

11 

11 

11 

19 

Measurement Points 

590 

462 

506 

539 

930 


Probe pressure measurements The flow field was measured with calibrated three-and 
five-hole probes. Data were accumulated in five measurement planes, each perpendicular to 
the duct centerline. The measurement plane locations along the duct centerline are presented 
in Table 113.1. These planes are also depicted by the shaded cross sections in Figure III.5. At 
each measurement plane the S-duct contained probe port holes, allowing access to the flow. 
The number of radial traverses and measurements made in each plane are given in Table 
in.l. Representative schematics depicting the location of probe port holes and measurement 
points are shown in Figures III.6 (Plane A) , III.7 (Planes B, C and D), and III. 8 (Plane E). 
Traversing intervals in the radial direction were approximately 0.254 cm for all measurement 
planes. In Planes A and E, traversing intervals near the surface were 0.0635 cm. Reported 
flow measurements were concentrated in only one symmetric half of the duct. 

A calibrated three-hole probe was used in Plane A, since a nearly one-directional velocity 
field existed there. The three-hole probe allowed total pressure and static pressure and two 
components of velocity to be measured. A yaw-nulling measurement method was utilized. 
The three-hole probe calibration method and data reduction technique are presented in 
Appendix A. 

A calibrated five-hole probe was used in Planes B, C, D and E, enabling total pressure, 
static pressure and three components of velocity to be measured. Both yaw-nulling and non- 
yaw-nulling measurement methods were employed. The five-hole probe calibration methods 
and data reduction procedures are presented in Appendix B. 

Surface static pressure measurements A total of 220 static pressure taps were 
located on the surface of the duct. The taps were constructed by locally drilling a 0.143 cm 
diameter hole normal to the duct surface. Each hole was plugged with an appropriate size 
aluminium tube having an inside diameter of 0.0813 cm. Epoxy placed on the outside of the 
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duct sealed the interface between duct and tube . The tube was sanded flush with the inside 
duct surface. A drawing representing a typical static tap is shown in Figure III.9. 

The static taps were distributed in the streamwise and circumferential directions. Three 
lines of taps ran in the streamwise direction. These streamwise lines were at constant angles 
of <f> = 10°, 90°, and 170°. Measurements at <f> = 0° and 180° were impossible because of 
the flanges there. The streamwise taps were spaced every s/d = 0.0873, beginning at s/d 
= 0.3492. A total of 53 taps were in each streamwise line. Four lines of taps ran in the 
circumferential direction. These circumferential taps were in Planes A, B, C and D. Plane 


z 

▲ 



Plane E 


Figure IH5 Half shell representation of the diffusing 
S-duct indicating the measurement plane locations 
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Figure III. 6 Measurement plane schematic for Plane A 
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Figure IIL7 Measurement plane schematic for Plane B, C and D 
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Figure III.8 Measurement plane schematic for Plane E 
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E was void of static taps. The circumferential static tap distribution is shown in Figure III. 6 
for Plane A and in Figure 113.7 for Planes B, C and D. 

Pressure measurement system A PSI Model 780B/T measurement system [26] was 
used to acquire all pressure data. The PSI system was a fully integrated measurement 
instrument consisting of electronically scanned pressure sensors and a microcomputer based 
data acquisition system. Individual pressure transducers provided high data acquisition rates 
for multiple pressure measurements. The accuracy of the system was maintained by frequent 
on-line calibrations of all transducers, performed transparently to the user at 20 minute 
intervals. After acquisition, information was carried by Escort, a data routing network, to a 
storage area for later post-processing on Sun SPARCstations. 

The PSI system contained four major components: the system controller, the data 
acquisition and control unit (DACU), the pressure calibration unit (PCU), and the sensor 


0.032 inch i.d. aluminium tube 
positioned normal to the surface 


duct surface 



epoxy 


Figure III.9 Static tap mounting procedure schematic 
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modules. The system controller interfaced the user to the DACU and directed the data flow 
within the system. The DACU managed the data acquisition functions and was controlled 
by an eight bit microprocessor. The DACU also interfaced with the PCU and reduced the 
transducer calibration data. The PCU applied the three point calibration to all transducers. It 
consisted of pneumatic valves and a highly accurate Digiquartz 103 KPa absolute pressure 
transducer (accurate to within + 0.02% full span). The Digiquartz transducer was used to 
calibrate the individual measurement transducers. Sensor modules were simply the link 
between the physical pressures and the electronic signals. Each module contained 32 
individual sub-atmospheric differential pressure transducers having a full scale range of 
34.5 KPa referenced to atmospheric pressure. The inaccuracy of an individual transducer 
was less than ± 0.10% of the full span. 
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CHAPTER IV 

COMPUTATIONAL ANALYSIS 

A second goal of this study was to numerically approximate the flow through the test 
S-duct. For the objectives of this thesis, a parabolized Navier-Stokes algorithm (PNS) was 
considered appropriate. In this chapter the formulation of governing equations, numerical 
procedure, grid generation and boundary conditions are described. 


Formulation of Equations 

Parabolized Navier-Stokes equations represent a mathematical model of three- 
dimensional viscous flows. PNS equations are intermediate in complexity between potential 
flow equations and full Navier-Stokes equations. PNS algorithms predict three-dimensional 
flow fields by solving a sequence of two-dimensional cross-stream governing equations 
while marching in the primary flow direction. This simplified solution procedure decreases 
computational time and storage requirements. 

For subsonic flows, the PNS equations are classified as partially-parabolized Navier- 
Stokes (P-PNS) equations. P-PNS equations are applicable for internal flows in which 
a predominant flow direction exists. P-PNS equations are derived from the steady Navier- 
Stokes equations by neglecting the effect of momentum diffusion in the streamwise direction. 
This approximation represents an attempt to make use of the experimentally observed fact 
that small disturbances at a given point, in a high Reynolds number ducted flow, are not 
transmitted very far up stream by viscous diffusion from that point. P-PNS equations also 
require that an initial approximation to the static pressure field is known a priori. This 
prerequisite incorporates elliptical influences, associated with subsonic flows, and static 
pressure gradients, associated with streamline curvature, into the solution. The initial static 
pressure distribution can be obtained from a solution to the velocity potential equation 
(Equation IV. 1) for a given duct geometry. Here, the local speed of sound, a, in Equation 
IV. 1 is given in Equation IV.2, where a 0 is constant. 
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P-PNS equations are valid in both the inviscid and viscous portions of the flow field, 
with the interaction between these regions automatically taken into account. One important 
restriction of the model is that the velocity component in the primary flow direction must 
be greater than zero. Hence, no flow reversal in the streamwise direction is permitted. 
No restrictions are placed on cross-flow velocities. Mathematically, partially-parabolized 
Navier-Stokes models can be represented by Equations IV.3— IV.8 for a Cartesian system 
where the x-component is in the direction of primary flow. 
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P = fi(p,T), P = f2(T), k = f 3 (T) (IV.9) 

The continuity equation, given in Equation IV. 3, is applicable at all points in the flow 
field. The momentum equations are represented by Equations IV.5— IV.7. Here, the static 
pressure at a point is the sum of, P, p m and p. The pressure, P, is the initial static pressure 
obtained from the potential flow solution. The components of its gradient are treated as 
source terms during the solution process. The mean streamwise pressure, p m , (Equation 
IV.5) is assumed to vary only in the streamwise direction and is determined with the aid 
of a global mass flow constraint (Equation IV.4). In contrast, the cross-stream pressure, 
p, in Equations IV.6 and IV.7, is permitted to vary across the channel. This uncoupling 
procedure assumes that small pressure variations across the channel have negligible effects 
on the streamwise momentum equation. The energy equation, given in Equation IV. 3, can 
be solved uncoupled from the momentum equations stated above for a constant stagnation 
enthalpy. Finally, state equations are needed to relate thermodynamic variables and transport 
properties as shown in Equation IV.9. 

Numerical Procedure 

The original source code, PEPSIG, used in this thesis was developed by Scientific 
Research Associates under NASA Lewis Research Center sponsorship. Modifications to 
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the code have been incorporated by Levy et al. [27], Briley and McDonald [28] and 
Anderson [29]. Further alterations to the code were performed by Wellborn [30], enabling 
the algorithm to run efficiently on UNIX workstations. Solutions to the PNS equations, 
by PEPSIG, were obtained with a two step process. A potential flow solution was first 
computed to approximate the static pressure field. Viscous calculations were performed next 
by marching once through in the primary flow direction. 

The potential flow numerical procedure involved solving the three-dimensional velocity 
potential equation, given in Equation IV. 1, including all cross derivatives on body 
fitted coordinates. A general purpose three-dimensional scalar finite-difference alternating- 
directional implicit (ADI) algorithm was used to solve Equation IV. 1. The three-dimensional 
approximate factorization technique of Douglas-Gunn was used to generate the ADI scheme 
from a basic Crank-Nicolson spatial differencing. The Douglas-Gunn scheme separated 
the three-dimensional operator into one-dimensional components and therefore split each 
iteration level into three steps, each step involving the implicit operations originating from 
a single coordinate. A simple tri-diagonal system resulted for each of the three steps. The 
potential function was then solved for iteratively until convergence. 

The viscous flow numerical procedure was based on the decomposition of the velocity 
field into primary and a secondary velocities. The primary velocity was determined from the 
streamwise momentum equation. The streamwise static pressure gradients were approximated 
from the inviscid flow solution. The one-dimensional, streamwise static pressure gradient 
correction was determined by ensuring global conservation of mass. These pressure gradients 
were treated as source terms in the streamwise momentum equation. The two-dimensional 
approximate factorization technique of Douglas-Gunn was used to create an ADI scheme 
which calculated the new primary velocity field. The secondary velocity was determined from 
a vector potential equation and a vorticity transport equation. A coupled ADI procedure was 
used to update these velocities. Corrections to these updated velocities were then computed. 
These correction velocities were assumed to be irrotational and ensured the difference form 
of the continuity equation was satisfied. Static pressures were then updated with the cross- 
stream momentum equations. FLARE approximations were used in regions of the flow 
field where the primary velocity was in the negative streamwise direction. The turbulent 
boundary layer was approximated by an algebraic turbulence model. The eddy mixing length 
was determined from an empirical relationship developed by McDonald and Camarata [31] 
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for equilibrium turbulent boundary layers. The equation of state for a perfect gas and a 
known constant stagnation energy closed the system of equations. 

The basic algorithm for PEPSIG is presented below. The input parameters are 
documented in Reference [32] and therefore not listed here. All governing equations were 
solved on a body fitted computational coordinate system. 

Inviscid Calculation 

• Generate a three dimensional grid 

• Solve velocity potential equation to obtain a static pressure field 

• Store pressure field 

Viscous Calculation 

• Generate a two-dimensional cross-stream grid 

• Solve primary momentum equation to update primary velocity 

• Solve vector potential and vorticity equations to update cross-channel velocities 

• Solve scalar potential equation to correct the cross-channel velocities 

• Solve Poisson pressure equation to correct the cross-stream pressures 

• Solve energy equation (if needed) 

• Step to next plane and start viscous calculation over again 

Grid Distributions 

Two computational meshes existed for the numerical solution. The grid distributions 
used contained the maximum number of grid points allowed by the limited size of the 
computer memory. Both meshes modeled the experimental S-duct and were made with an 
algebraic grid generator. 

The potential solution utilized 50 x 19 x 10 grid points in the stream wise, circumferential 
and radial directions. An O-grid was implemented in the cross-stream plane. The O-grid 
was equally spaced in radial and circumferential directions. The centerline coordinates of the 
experimental S-duct were defined with a tenth degree polynomial which was a least-squares 
approximation of Equations in.l, III.2 and III.3. The variation of radius in the stream wise 
direction was given by Equation III.4. The cross-sectional and streamwise distributions for 
the potential grid are shown in Figure IV. 1. 


a) Cross-sectional view 



b) Streamwise view 

Figure IV. 1 Computational mesh distribution for the potential solution 
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Symmetry 



a) Cross-sectional view 



Figure IV.2 Computational mesh distribution for the viscous solution 
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The viscous solution utilized 77 x 49 x 49 grid points in the streamwise, circumferential 
and radial directions. Again, an O-grid was implemented in the cross-stream plane. However, 
stretching was used in both the circumferential and radial directions. In the radial direction, 
clustering near the duct surface was used to resolve the boundary layer. In the circumferential 
direction, clustering near $ = 180° was used to help enhance the calculations in the separated 
flow region. The variation of radius in the streamwise direction is given by Equation III.4. 
The streamwise positions were defined by the same polynomial used for the potential solution. 
The cross-sectional mesh and a representative streamwise distribution for the viscous grid 
are shown in Figure IV.2. 


Boundary Conditions 

Two sets of boundary conditions were used for the numerical solutions. The potential 
solution met a flow tangency condition on the duct surface. Zero gradients of the 4> function 
along the plane of symmetry were also satisfied. To achieve zero gradients, a reflective 
boundary at the plane of symmetry was assumed. The centerline was considered a pole 
boundary. A uniform flow distribution existed at both the inlet and exit. The downstream 
potential distribution, <3>, was specified. These conditions are indicated in Figure IV. 1. 

For the viscous solution, no slip on the duct wall and zero gradients along the plane of 
symmetry were used. Again, a reflective boundary at the plane of symmetry was assumed 
and the centerline was considered a pole boundary. The centerline Mach number, Reynolds 
number, static temperature, static pressure and the initial boundary layer thickness at the first 
marching station were specified. These conditions are indicated in Figure IV.2. 
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CHAPTER V 

INLET FLOW CONDITIONS 


Test inlet flow conditions were obtained from a survey of the flow field one radius 
upstream of the first bend (Plane A in Figure III.5). All reported tests were conducted with 
an inlet centerline Mach number of 0.6. The Reynolds number, based on the inlet diameter 
and centerline velocity, was 2.6 x 10 6 . Hot wire measurements were previously gathered by 
Reichert [33] to ascertain the turbulence intensity at the centerline in Plane A. Measurements 
were made for Mach numbers ranging from 0.05 to 0.5 and for five overheat ratios. A 
50Khz data acquisition sampling rate was used. The centerline turbulence intensity in Plane 
A for Reichert’s study [33] was nominally 0.65%, which was presumed to represent the inlet 
centerline turbulence intensity in this thesis. 

A thin turbulent inlet boundary layer existed in Plane A. The boundary layer thickness, 
6, was defined as being from the wall to where 95% of the free stream velocity was 
achieved. The displacement thickness, <5i, momentum thickness, <$ 2 , and shape factor, H, 
used to further quantify the inlet flow conditions, are defined by Equations V.l— V.3. These 
boundary layer parameters were calculated by numerically integrating the survey data. The 
inlet flow conditions are listed in Table V.l. 


61 = 



6 2 = 



(V.l) 


(V.2) 



(V.3) 


A representative inlet velocity profile, plotted in nondimensional law-of-the-wall 
coordinates, is shown in Figure V.L The friction velocity, used to define law-of-the-wall 
coordinates, was not measured but instead was chosen to provide the best approximation 
of the linear profile region of the test data to the law-of-the-wall logarithmic function. 
Comparisons indicate little deviation from a conventional turbulent boundary layer. 
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Table V.l Experimental and computational flow conditions in Plane A 


Inlet Parameter 

Experimental 

Computational 

M d 

0.60 

0.60 

Re d 

2,600,000 

2,600,000 

(6/r i) x 100 

6.95 

4.45 

(<5i/n) x 100 

1.46 

0.71 

(S 2 /r i) x 100 

1.06 

0.53 

H 

1.38 

1.35 


The numerical solution started two inlet diameters upstream of Plane A. Therefore, 
computational input parameters corresponded to flow variables at that position while test 
flow conditions pertained to flow variables in Plane A. This posed a problem when defining 
the computational input parameters. In practice it was very difficult to chose the input 
parameters which gave an exact match of all computational flow conditions in Plane A to 
the test flow conditions. The numerical and test flow values of Mach and Reynolds numbers 
in Plane A were equal. However, exact matching of computational and experimental flow 
boundary layer flow parameters in Plane A was unfortunately not attained. The computational 
flow conditions for Plane A are listed in Table V.l. 



Figure V.l Experimental law-of-the-wall velocity profile measured at Plane A 
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CHAPTER VI 

RESULTS AND DISCUSSION 


All data presented in this chapter are in nondimensional form. Aerodynamic results 
represent mean values of measurements. Pressures are presented as total and static pressure 
coefficients, which are defined by Equations VI. 1 and VI.2. The pressures p 0 and p represent 
local values of total and static pressures. Inlet centerline conditions define the reference 
states p 0jC i and p c /. Three dimensional velocity components were converted to local Mach 
vectors and normalized by the inlet centerline Mach number, as shown in Equation VI.3. 
Finally, streamwise vorticity, given in Equation VI.4, was normalized by the ratio of inlet 
centerline velocity to inlet duct diameter. 


Po - Pci 
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(VI. 1) 
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For this section the nondimensional pressure coefficients, normalized Mach vector and 
normalized streamwise vorticity component will simply be referred to as total pressure, 
static pressure, velocity and vorticity. The velocity vector will be presented as normal and 
transverse components. The normal component is perpendicular to the measurement plane 
while the transverse component is parallel to the measurement plane. Transverse velocity 
components, originally on a polar grid, were interpolated onto a Cartesian mesh [33], in 
order to enable easy visualization and comparison of the cross-stream flow fields. Since data 
were measured and computed in only one half of the symmetrical duct flow, all cross-stream 
aerodynamic results were mirror imaged to enhance presentation. 
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Surface Streakline Patterns 

Major test flow field characteristics near the duct wall were visualized with streaklines of 
fluorescent oil on the duct surface. Three conclusions were made from the streakline patterns. 
First, the flow was indeed symmetric. This was ascertained by applying a different color of 
dye to the surface of each symmetric duct half. The dye from each half remained separate 
and identical streakline patterns were present on each symmetric half of the duct. Symmetry 
is illustrated in Figure VI.2. Here, one symmetric half of the duct is shown connected to the 
exit extension. The flow is from left to right. Obviously, red dye from the left half of the 
duct does not mix with the green dye from the right half of the duct. This is true everywhere 
except along the vertical split line where there exists a thin line of yellow dye (mixed red and 
green dyes). Second, a large region of three-dimensional separated flow existed, as strikingly 
indicated in Figures VI.3 and VI.4. The entire duct section is shown in Figure VI.3, while 
Figure VI.4 is a close up of the separated flow region. The free stream flow is from left 
to right and each photo shows only one symmetric half of the observed streakline pattern. 
The separated flow region consisted of two saddle points, occurring on the duct split line 
(</> = 180°), and two spiral nodes, lying in each symmetric duct half. The entire separated 
region was located on the lower duct wall. The onset of separation (upstream saddle point) 
was located at s/d = 2.02, while reattachment (downstream saddle point) occurred at s/d = 
4.13. Three-dimensional separated flow terminology is explained in detail in Chapter VII. 
Third, boundary layer cross flows were present (Figure VI.3). In the first bend, streaklines 
were driven toward the lower surface ( <f> = 180°). In the second bend, streaklines near 
</> = 90° diverged. The upper streaklines converged toward the top of the duct (</> = 0°), 
while the lower streaklines were directed toward the bottom of the duct (d> = 180°). 

A qualitative assessment of the flow in the plane of symmetry ( x-z plane) was acquired 
by temporarily introducing a splitter plate between the two symmetric duct halves. The plate 
was installed for only one test to produce streakline patterns on the plate. The S-duct surface 
results given above were obtained without the splitter plate in place. The splitter plate results 
are assumed to be qualitatively correct, even though the presence of the plate introduced 
additional shear into the flow. The streakline patterns produced on the splitter plate are 
presented in Figure VI.5. Free stream flow is from left to right. The location of separation 
and reattachment are identifiable. The streaklines follow the duct curvature in the first bend. 
Blockage caused by flow separation forced the streaklines to deviate from the curvature of 
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Figure VI.2 Surface oil flow patterns at the duct exit showing flow symmetry 



Figure VI.3 Surface oil flow patterns on one symmetric half of the duct surface 
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Figure VI.4 Close up of the surface oil flow patterns developed in the separation region 



Figure VI.5 Surface oil flow patterns on the centerline splitter plate 
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the second bend. The reversed flow region began thin and grew until reattachment. The 
extensive bending of the streaklines at the exit of the duct was caused by viscous effects 
associated with the splitter plate. Cross-stream pressure gradients, at the duct exit, helped 
drive the boundary layer fluid on the splitter plate toward the top of the duct. This was most 
noticeable at the exit where the splitter plate boundary layer was thickest. 

Surface Static Pressures 

Symbols and lines in Figure VI.6 represent the experimental and computational surface 
static pressure variations with axial distance for three circumferential locations. The region 
of streamwise separated flow, deduced from flow visualization, and the corresponding 
computational separated flow region are also shown. 

For the experimental measurements, the effects of streamline curvature and diffusion are 
clearly indicated within the first bend by the pressure difference between the <j> = 10° and 
the <f> = 170° data and the overall pressure rise. The influence of flow separation is shown 
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Figure VI.6 Axial distributions of surface static pressure for three circumferential positions 
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by the constant values of the <f> = 90° and 170° static pressures between 2.0 < s/d < 3.2. 
The proximity of these two curves, between 2.0 < s/d < 3.2, indicate only small changes in 
flow speed and/or direction in the lower half of the duct. Blockage caused by the separated 
flow also increased the favorable pressure gradient at <f> = 10°. Static pressures at the three 
circumferential positions are nearly equal at s/d = 2.6. This suggests that a uniform static 
pressure distribution exits throughout the cross-stream plane there. Static pressure rose again 
for s/d > 3.2, even though the flow remained separated beyond this point, because the duct 
continued to diffuse the flow. After reattachment the pressure distributions converged to 
Cp = 0.466 far downstream (. s/d « 9.0). 

Computational results show the same distribution trends mentioned above for the static 
pressure measurements. Upstream of flow separation, streamline curvature and diffusion 
effects were evident. The computed static pressures at the three positions (10°, 90°, 170°) 
were nearly equal at s/d = 2.6, as with the test measurements. Downstream of the separated 
flow region, the pressure values converged to a single value, Cp = 0.517. The major 
discrepancies between computed and experimental results occurred within the second bend 
of the duct. The computed region of separated flow, 3.15 < s/d < 4.89, was located further 
downstream than the experimental region. Also, the computed static pressure values at 
<j> = 90° and 170° did not coincide in this region as did the measured values. Numerical 
results of the static pressure distribution in this S-duct, reported by Harloff et al. [19], 
have been obtained using both algebraic and k-t turbulence models with a FNS algorithm. 
The FNS results using both turbulence models within the second bend were similar to the 
P-PNS calculations. Therefore, errors caused by simplifications to the streamwise static 
pressure within the P-PNS equations are believed to be negligible when compared to the 
errors caused by the choice of turbulence models. 

The circumferential distributions of surface static pressure in Planes A through D are 
presented in Figure VI.7. Symbols represent the experimental data. Lines indicate the 
computational results. The experimental and computational pressures for Plane A were 
nearly equal, both circumferentially and to each other. Essentially no influence of the duct 
curvature on the flow was detected at this upstream plane. The static pressure data in Figure 
VI.7 for Planes B, C and D all reflect the influence of streamline curvature. The pressures 
measured in Plane B were largest at <t> = 10° and continually decreased until <j> = 170°. The 
values were negative for <f> > 110°. Separation did not affect the flow field in Plane B. The 
distribution of computational and experimental static pressures at Plane B were similar. In 
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Planes C and D the cross-stream pressure distribution reversed from what it was in Plane 
B. Test measurements show that the separated flow reduced the static pressures at large 
values of <j> in Planes C and D. The peak pressure location was shifted to a larger value 
of <f> between Planes C and D. This trend indicates a reduction in size of the region of 
separated flow from Plane C to D. The computed pressure distributions, in Planes C and 
D, agreed in shape but were overall higher in level than the experimental pressures. This 
was attributed to two factors; first, the delayed prediction of separated flow and second, less 
computed blockage than actually occurred. 

Cross-Stream Aerodynamic Data in Plane A 

The flow field in Plane A corresponded to developing pipe flow. The results in Figure 
V.l and Table V.l verify the that boundary layer was fully turbulent. Total pressure contours 
for Plane A are presented in Figure VI.8. Normal velocity contours are displayed in Figure 



Figure VI.7 Circumferential distributions of surface static pressures in Planes A, B, C and D 
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a) Experiment 



Figure VI. 8 Total pressure distributions in Plane A 
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a) Experiment 



b) Computation 

Plane A s / d = -0.50 
E 

Figure VI.9 Normal velocity distributions in Plane A 
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VI.9. Expanded views of these distributions near <j> = 90° are included for comparison. 
Both the normal velocity and total pressure distributions (experimental and computational) 
show no circumferential variance. The measured and computed cross-stream static pressure 
distributions were nearly uniform, and are therefore not shown. This uniformity is consistent 
with the previously reported surface static pressures (Figures VI.6 and VI.7). The measured 
and computed cross flows were negligible, so transverse velocities are not presented. Vorticity 
was present only in planes normal to the streamwise direction. 

The only discrepancy between the experimental and computational results is with the 
boundary layer thickness. This difference was explained in the Inlet Flow Conditions Chapter. 
The computational boundary layer was thinner than the experimental boundary layer. Since 
centerline and not bulk conditions were used to normalize aerodynamic data, this discrepancy 
does not affect the legitimacy of downstream comparisons. 

Cross-Stream Aerodynamic Data in Plane B 

Initially the S-duct deflected the incoming flow downward. This caused the flow field to 
deviate significantly from traditional developing pipe flow by Plane B. Cross-stream static 
pressure gradients developed as the flow progressed downstream. These gradients were 
a direct result of the core flow adjusting to duct geometry and the resulting streamline 
curvature. The measured and computed static pressure distributions in Plane B (Figure 
VI. 10) indicate that the maximum static pressure existed at the top of the duct (<fr = 0°) and 
the minimum was at the bottom of the duct ( <j> = 180°). These distributions are consistent 
with the surface static pressure data in Figures VI.6 and VI.7. 

Total pressure contours are shown in Figure VI.ll. Numerical and experimental results 
are nearly identical. There was little change in the total pressure distribution from Plane A 
to Plane B (Figures VI.8 and VI.ll) except for a slight thickening of the boundary layer 
with downstream distance. The boundary layer thickness in Plane B did vary slightly with 
circumferential position. The boundary layer thickness was greatest at <f> = 0° and least at 
<f> = 180°. This trend can be attributed to the streamwise pressure gradients upstream of 
Plane B which accelerated the flow near the bottom of the duct and decelerated the flow 
near the top of the duct. 

The normal velocity distributions in Plane B (Figure VI. 12) reflect the influence of 
streamline curvature. To reiterate, at the top of the duct, where the cross-stream pressure 
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a) Experiment 



Figure VI. 10 Static pressure distributions in Plane B 
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a) Experiment 



Figure VI. 11 Total pressure distributions in Plane B 
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a) Experiment 



Figure VI. 12 Normal velocity distributions in Plane B 
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a) Experiment 



Figure VI. 13 Transverse velocity vectors in Plane B 
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a) Experiment 



Figure VI. 14 Stream wise vorticity distributions in Plane B 
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was largest, an adverse streamwise pressure gradient decelerated the flow. At the bottom of 
the duct, where the minimum cross-stream pressure existed, a favorable streamwise pressure 
gradient accelerated the flow. This acceleration can be seen in Figure VI. 12, where the local 
Mach number near the bottom of the duct exceeded the inlet Mach number. 

Experimental and computational transverse velocity components in Figure VI. 13 show 
that cross flows occurred only near the duct surface. These cross flows were caused by the 
pressure gradient turning the lower momentum boundary layer fluid toward the bottom of the 
duct. This was also detected in flow visualization, as seen in Figure VI.3. It is interesting 
to note that even though the pressure-driven cross flows near the duct wall were occurring 
in Plane B, a large region of low momentum fluid was not present near the bottom of the 
duct (Figure VI. 11). This changed by Plane C. 

The vorticity plots in Figure VI. 14 confirm that transverse velocities in Plane B were 
confined to the wall boundary layer region only. Here the solid lines indicate positive 
vorticity, while dashed lines represent negative vorticity. The largest vorticity magnitude 
and hence the most amount of turning occurred at the sides of the duct. This was also 
illustrated in the transverse velocity plot (Figure VI. 13). Excellent agreement in contour 
shapes and magnitudes of streamwise vorticity existed between test and computational data. 

Cross-Stream Aerodynamic Data in Plane C 

The location of Plane C (s/d = 2.97) was 0.35 diameters downstream of the middle of the 
S-duct. This plane was in the midst of the actual separated flow region. Differences between 
the computed and measured aerodynamic data arose since the computed separated flow region 
occurred further downstream (s/d = 3.14). The contours of static pressure, total pressure 
and normal velocity for Plane C are shown in Figures VI. 15, VI.16 and VI.17. Transverse 
velocity components and streamwise vorticity are depicted in Figures VI. 18 and VI. 19. 

Near the middle of the S-duct the static pressure distribution should become nominally 
flat as the centerline curve undergoes an inflection. In the second bend of the duct the 
orientation of the cross-stream static pressure distribution should reverse, so that the lowest 
static pressure is near the top of the duct and the highest static pressure is near the bottom. 
The static pressure data for Plane C in Figure VI. 15 confirm this reversal. The measured 
static pressure data in Plane C also revealed a local region of nearly constant static pressure 
in the lower half of the duct. This trend substantiates a previous conclusion drawn from 



a) Experiment 



Figure VI. 15 Static pressure distributions in Plane C 
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a) Experiment 



Figure VI. 16 Total pressure distributions in Plane C 
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a) Experiment 



Figure VI. 17 Normal velocity distributions in Plane C 
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Figure VI. 18 Transverse velocity vectors in Plane C 
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a) Experiment 



Figure VI. 19 Streamwise vorticity distributions in Plane C 
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the surface pressure data, namely, the existence of a region of only small changes in flow 
speed and/or direction. It is also consistent with the separated flow there. This same trend 
was not demonstrated by the computational results because of the delayed prediction of 
flow separation. 

Strong, static-pressure-driven boundary layer cross flows existed near the lower duct 
surface in Plane C (Figure VI. 18). These cross flows continually moved boundary layer 
fluid into the low velocity region. The cross flows suggest that two vortical structures were 
developing. These vortices were not present in Plane B. A significant amount of upward 
deflected flow at the centerline in Plane C is also shown in the measured data. This confirms 
the results obtained from the surface flow visualization on the temporary centerline splitter 
plane. This deflection was a direct consequence of two factors, blockage created by flow 
separation and cross flows driven by pressure gradients. Computationally, the smaller amount 
of upward centerline velocity deflection was a result of pressure-driven cross flows only. 

Total pressures (Figure VI. 16) and normal velocity components (Figure VI. 17) show 
the accumulation of low momentum fluid in the bottom portion of the duct. The measured 
accumulation was caused by the reversed flow, the adverse streamwise pressure gradient and 
the convection of boundary layer fluid toward the lower surface. The same accumulation for 
the computational results was caused only by the adverse streamwise pressure gradient and 
the cross passage convection of boundary layer fluid. 

The streamwise vorticity data, shown in Figure VI. 14, again confirm the above mentioned 
trends. Here, the position of the largest vorticity magnitude moved slightly lower but was 
still at the sides of the duct. Overall agreement between test and computational data of the 
magnitudes of vorticity was surprisingly good, except near the bottom of the duct. This was 
where actual separation was present. 

Cross-stream Aerodynamic Data in Plane D 

The location of Plane D was at s/d = 4.01. This position was close to the actual 
reattachment point and in the middle of the computed separated flow region. Again, 
differences between the computed and measured aerodynamic data were expected and are 
present. For Plane D, the contours of static pressure, total pressure and normal velocity are 
shown in Figures VI.20, VI.21 and VI.22. Transverse velocity components are depicted in 
Figure VI.23. Streamwise vorticity distributions are plotted in Figure VI.24. 
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Figure VI.20 Static pressure distributions in Plane D 
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a) Experiment 



Figure VI.21 Total pressure distributions in Plane D 
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a) Experiment 



Figure VI.22 Normal velocity distributions in Plane D 
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a) Experiment 





Figure VI.23 Transverse velocity vectors in Plane D 


57 



a) Experiment 



Figure VI.24 Streamwise vorticity distributions in Plane D 
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The static pressure data in Plane D (Figure VI.20) are similar to static pressure data 
in Plane C. The measured distribution still indicated a substantial region of nearly constant 
static pressure in the lower half of the duct (Cp ~ 0.35). This region was also present 
in Plane C. The computational results suggest a similar region, although at a much higher 
pressure value ( Cp « 0.475). 

Total pressures (Figure VI.21) and normal velocity components (Figure VI.22) suggest 
further growth of the low total pressure/low velocity region in the bottom half of the duct. 
Measurements indicate that enough low momentum fluid was carried away from the duct 
wall to form a region of near zero velocity far from the surface. Computationally this region 
was still evolving. 

Strong pressure-driven cross flows still existed near the lower duct surface in Plane 
D (Figure VI.23). These cross flows continually drove boundary layer fluid toward the 
low velocity region. A strong vortical structure developed. For the experiment, the core 
flow returned to the nominal streamwise direction by the time it reached Plane D, which 
is contrary to the splitter plate streakline patterns. This inconsistency was caused by the 
pressure gradients which drove the boundary layer fluid on the splitter plate toward the top 
of the duct, as previously explained. 

Streamwise vorticity data (Figure VI.24) revealed that the boundary layer fluid at the 
top of the duct began to turn upward. This was slightly noticeable in the transverse velocity 
plots (Figure VI.23). This reversal in cross-stream flow direction was caused by the change 
in static pressure distribution from Plane B to Plane D. It is interesting to note that even 
though Plane D is well into the second bend, no large vortical structure existed at the top 
of the duct in Plane D. 

Cross-stream Aerodynamic Data in Plane E 

By Plane E the free stream flow returned nominally to the ^-direction. Cross-stream 
static pressure gradients were nearly eliminated and are therefore not illustrated. The total 
pressure and normal velocity distributions are represented by contours in Figures VI.25 and 
VI.26. Transverse velocity components and streamwise vorticity for Plane E are depicted 
in Figures VI.27 and VI.28. 

Even though cross-stream static pressure gradients did not exist in Plane E, upstream of 
Plane E, pressure gradients were strong enough to fully reverse the direction of the boundary 
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a) Experiment 



Figure VI.25 Total pressure distributions in Plane E 
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Figure VI.26 Normal velocity distributions in Plane E 
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a) Experiment 



E 


Figure VI.27 Transverse velocity vectors in Plane E 
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a) Experiment 



Figure VI.28 Streamwise vorticity distributions in Plane E 
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layer fluid flow on the upper surface of the duct, as seen in Figures VI.27 andVI.28 for 
both measured and computed results. This was also detected with flow visualization. This 
phenomenon had just begun in Plane D. A pair of counter-rotating vortices, located in the 
lower half of the duct, evolved. These vortices continually convected the low momentum 
fluid of the boundary layer toward the center of the duct. The positions of the computed 
counter-rotating vortices were not in the same location as the measured vortices. This 
discrepancy between computed and experimental results was noted also by Harloff et al. 
[19]. Vorticity values (Figure VI.28) also indicated that the measured vortical structure was 
much weaker than the computed vortical structure. 

Low momentum fluid convection by the vortices degraded both the uniformity and 
magnitude of the exit total pressure distribution. The measured region of low momentum 
fluid in Plane E (Figure VI.25 and VI.26) extended above the duct center line. This trend had 
been observed before by previous investigators [5,9], The inability of numerical algorithms 
to capture this flow phenomenon (as in the P-PNS results) has also been documented by 
several researchers [14,18,19]. 
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CHAPTER VII 


FURTHER INSIGHTS INTO SEPARATED FLOW 


Flow separation is usually related to the occurrence of reversed flow near a solid surface. 
Often, reversed flow simply refers to any flow whose direction is opposite to the projection 
of the free stream direction onto the surface. Therefore, in a two-dimensional steady flow, 
the phenomenon of flow separation may be clearly defined. Unfortunately, two-dimensional 
flow is very rare. For three-dimensional flows, the specification of separation in terms of 
reversed flow may be inadequate and/or inappropriate. This chapter attempts to describe and 
classify the region of separated flow within the S-duct by using well defined terminology. 
This terminology has its origin in notions based on topological structures and streamsurface 
bifurcation. 

Legendre proposed that a pattern of streamlines near a surface may be considered 
as trajectories having properties consistent with those of a vector field [34], A principle 
property of a vector field states that through any nonsingular point there must pass one and 
only one trajectory. Therefore if streamline trajectories near a surface constitute a vector 
field, elementary singular points of this field can be categorized mathematically. Departing 
from the notion of nearwall streamlines, whose velocities vanish to zero at the wall, Lighthill 
[35] connected skin-friction lines on a body surface to this continuous vector field. In either 
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Figure VII. 1 Illustrations of elementary singular points 
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case, it is common to assume that oil-streakline patterns which develop on a solid surface 
during an experiment closely represent both skin-friction lines and nearwall streamlines. 

Singular points in a pattern of skin-friction lines or oil-surface streaklines occur at points 
on the surface where skin-friction becomes identically zero. Two types of singularities can 
occur: nodes and saddle points. Nodes are further classified as nodal points or foci. For this 
thesis, foci are simply termed spiral nodes. 

A point common to an infinite number of skin friction lines is termed a nodal point. All 
but one of these skin-friction lines are tangent to a specific single line emanating from the 
singular point. An exception to this case is when every skin-friction line leaving or entering 
a singular point has a distinct tangent. This type of nodal point can be produced by flow 
over a thin axisymmetric cylindrical body. A spiral node has no common tangent line and 
therefore differs from a nodal point. In a spiral node, an infinite number of skin-friction 
lines spiral around the singular point. If the skin-friction lines spiral away from the point 
it is called a spiral node of attachment. If the skin-friction lines spiral into the point it is 
called a spiral node of separation. A saddle point differs from any node in that only two 
distinct skin-friction lines pass directly through the singular point. All other skin-friction 
lines around a saddle point miss the singular point. In this respect, a saddle point commonly 
acts as barrier between two adjacent sets of skin- friction lines. Representative illustrations 
of these singular points are presented in Figure VII. 1. 

With the notions of streamsurface bifurcation theory, Homung and Perry [36] have 
added to the vocabulary of well defined terms for the description of three-dimensional steady 
separated flows. Many steady three-dimensional flows can be described by streamsurfaces. 
Often in flows fields there exist special streamsurfaces which appear to bifurcate (split apart or 
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Figure VII.2 Streamsurface bifurcations on a solid wall 
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combine together). There are two ways a streamsurface can bifurcate. A positive bifurcation 
refers to a stagnation streamsurface which divides along a stagnation line. For a positive 
bifurcation, a single streamline within the stagnation streamsurface appears to bifurcate into 
two different streamlines at the positive bifurcation line. Positive bifurcations usually are 
associated with downwashing counter-rotating vortices. A negative bifurcation refers to a 
streamsurface which emerges from a negative bifurcation line. In a negative bifurcation, 
two different streamlines near a flow boundary combine to form a single streamline within 
the streamsurface. Negative bifurcations usually are connected to uplifting counter-rotating 
vortices. Positive and negative bifurcations are illustrated in Figure VII.2. 

From the few notions of streamsurface bifurcation and topological singularities introduced 
above, the region of separated flow within the S-duct can now be adequately described. To 
begin, imagine cutting the surface of the test S-duct along the top of the vertical centerline 
plane ( <f> = 0°). Now imagine unfolding the duct onto a horizontal plane. The oil-streakline 
patterns developed during the S-duct surface flow visualization tests (Figures VI.2— VI.5) 
then would closely resemble Figure VII. 3. Here the flow is from left to right and both 
symmetric halves of the duct surface are shown. These original streakline patterns can be 
mapped onto a “skeleton” drawing, which illustrates the important topological singular points 


Figure VII.3 Unwrapped S-duct surface streakline patterns 
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and bifurcation lines. The skeleton drawing related to the S-duct surface streakline patterns 
is shown in Figure VII.4. Two saddle points, two spiral nodes, two positive bifurcation lines 
and four negative bifurcation lines are identified in the simplified skeleton drawing. 

The streakline patterns indicate that two negative bifurcation lines emanated from the 
upstream saddle point. These negative bifurcation lines ran in opposite cross-stream directions 
and were symmetric about the centerline. Each of these bifurcation lines were drawn into 
two symmetrically arranged spiral nodes of separation. A third negative bifurcation line 
developed from the downstream saddle point. A fourth negative bifurcation line was located 
at the top of the duct (along the line which was cut to unfold the S-duct). A pair of positive 
bifurcation lines flanked the two spiral nodes on the downstream side of the pattern. 

This pattern of streaklines has been previously named an "owl face" separation [37]. 
Streakline patterns similar to this one have been observed in regions where the flow is under 
the strong influence of two uplifting vortices. These strong uplifting vortices were shown 
to exist in the S-duct. Owl face separations are further classified as first or second kinds. 
Owl face separations of the second kind have another spiral node within the flow field 
located on the symmetric plane . An owl face of the first kind is void of this spiral node 



Figure VII.4 Skeleton schematic of the surface streakline patterns 
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on the symmetric plane. After observing to the flow visualization patterns developed on the 
surface of the temporary splitter plate, the separation region can now be further classified 
as an “owl face of the first kind”, since no visible spiral node was present on the centerline 
splitter plate. A perspective drawing representing an owl face separation of the first kind 
is shown in Figure VII.5. The two tubes in Figure VII.5 correspond to the stream surfaces 
coiled around the two vortices. This drawing corresponds well to the S-duct experimental 
transverse velocity measurements and surface flow visualization observations. 

In Figure VII.4 a single streakline has been drawn between the two saddle points without 
passing through a node on the way. Some researchers [34] consider this to be an impossible 
situation because it is unstable. These researchers proclaim that slight variations in any flow 
field conditions would alter this topological structure. However, Perry and Homung [37] 
ingeniously used these slight variations in flow field conditions to justify why this topological 
structure can exist. They explain that the pattern drawn in Figure VII.4 is extremely special, 
and the slightest asymmetry in the flow would alter the streakline pattern to the one shown 
in Figure VII.6. Here, no trajectories pass directly between two saddles. The slightest 
unsteadiness might even cause the pattern to oscillate between that of Figure VII.6 and 
its mirror image. Perry and Homung further point out, that from a practical viewpoint the 



Figure VII.5 Perspective view of an owl face separation of the first kind 
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patterns in Figures VII.4 and VII. 6 are not experimentally distinguishable since the streakline 
emanating from a saddle point is almost impossible to accurately identify. 

Upon describing the separated flow region within the S-duct, the question arises, “What 
portion of the flow field is truly part of the separated flow?” and/or, “Is separated flow even 
a proper term to use when dealing with three-dimensional flow ?”. These questions are not 
readily answered. It is easy to see though that referring to separated flow strictly as reverse 
flow is incorrect, for within the S-duct only a small region between the two spiral nodes 
involved reversed flow. However, most aerodynamicists would agree that the flow in the 
area surrounding the spiral nodes and even perhaps the flow near the negative bifurcation 
lines constitutes separated flow. Many researchers deem that a necessary condition for the 
occurrence of flow separation on a continuous surface is the convergence of oil-streakline 
patterns onto a particular streakline. This is perhaps the best current indicator of three- 
dimensional separated flow, and it is supported by the present S-duct aerodynamic and 
surface flow visualization data. 



Figure VII.6 Skeleton schematic illustrating asymmetric surface streakline patterns 
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CHAPTER Vffl 
CONCLUSIONS 

Compressible, subsonic flow through a diffusing S-duct was experimentally and com- 
putationally investigated. Benchmark aerodynamic measurements of the three dimensional 
velocity field and total pressures and static pressures were obtained in five cross-sectional 
planes. Surface static pressure and flow visualization data were also acquired. The S-duct 
flow field was computed using a subsonic parabolized Navier-Stokes code. Experimental 
and numerical results were compared. 

The experimental measurements show that a large three-dimensional region of separated 
flow occurred within the duct. A large portion of the duct cross-section was blocked by this 
separated flow region. This blockage forced the core velocity to increase while passing over 
the separated flow region. The increase in core velocity led to larger shear stress levels and 
therefore larger total pressure losses within the S-duct. The duct curvature induced strong 
pressure driven secondary flows, which evolved into two counter-rotating vortices. These 
vortices convected low-momentum boundary layer fluid into the lower portion of the duct and 
from there out toward the center of the duct, degrading both the uniformity and magnitude 
of the total pressure distribution. At the duct exit, a large region of low momentum fluid 
filled the lower half of the duct and extended above the centerline. 

The computational results are generally in agreement with the experimental measure- 
ments. Strong cross flows, total pressure nonuniformity and a separated flow region were 
all predicted. However, the separation region was calculated to occur further downstream 
than actually observed. The calculated cross flows did not convect as much low momentum 
fluid to the center of the duct as the actual cross flows did. These discrepancies, common 
to the use of both full and parabolized Navier-Stokes algorithms to simulate duct flows, are 
linked mainly to turbulence modeling and were not necessarily the fault of the simplified 
governing equations used in parabolized codes. In fact, parabolized codes are better suited 
for preliminary design purposes because of inherent reduced computational effort. 

These experimental and computational analyses of the flow field could ultimately be 
beneficial to S-duct designers. The numerical results demonstrate how approximate flow 
field trends can be predicted using a simplified parabolized model of the flow. The three- 
and five-hole probe calibration and data reduction procedures could enable an experimentalist 
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to obtain quick and accurate measurement data in other ducts. The collected experimental 
data are suitable for the validation of computational codes. Data for S-ducts having improved 
performance could be compared to these baseline data. 
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CHAPTER IX 

RECOMMENDATIONS FOR FURTHER RESEARCH 

Although the material reported in the thesis is complete, additional experimental and 
computational analyses of the flow through a diffusing S-duct would be beneficial to inlet 
designers. Further baseline experimental studies using hot wires to obtain Reynolds stresses 
and an estimate of velocity fluctuations are still needed. Hot film gages, located on the 
surface of the duct, may be able to better trace the separation footprint, including some 
unsteady aspects of this complex flow phenomenon. In addition, an ethylene trace gas study 
would help confirm fluid particle pathline conclusions. 

Investigations involving active or passive flow management devices, which control 
the extent and effects of flow separation, would be very insightful. A parametric study 
determining the performance benefits acquired when embedded counter-rotating vortex 
generators are installed in the duct is currently under progress. A study of the effects 
of embedded co-rotating vortex generators or boundary layer suction on the flow field to 
manage separation would also be useful. 

The advancement of military aircraft continually make the designer expand current inlet 
flow technology. Advances in canard control-surfaces bring to light the need for experimental 
data for S-duct flows with asymmetric inlet profiles. Simple vortex ingestion experiments 
could produce valuable insights to this design problem. A second type of S-duct, having 
square-to-circular transitioning cross sections, is often used in military aircraft. Studies, 
similar to the one described herein need to be done on this type of S-duct. 

The conclusions of this thesis indicate that much effort must be channeled toward the 
development of an accurate anisotropic turbulence model in order to accurately predict 
complicated flow fields, as in an S-duct. Unfortunately, this is not an easy task. Therefore, 
designers and researchers must use the tools available, keeping in mind the limitations and 
strengths involved. 
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CHAPTER X 
POSTSCRIPT 

The aerodynamic benchmark data collected in this study may be used for the validation 
of computational fluid dynamic codes. Further information concerning this study may be 
obtained from the authors. The data is available on diskette, tape or through the internet. 
Please send internet mail to fsbud@lerc.nasa.gov or fsbud@iastate.edu . 
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APPENDIX A 

THREE-HOLE PROBE CALIBRATION AND DATA REDUCTION TECHNIQUES 

This appendix contains supplementary information on the calibration and data reduction 
of a three-hole probe. The information covering five-hole probes can be found in Appendix 
B. Both the three- and five-hole probe calibration techniques used in this thesis closely 
follow the methods developed by Reichert [33]. The reader is referred to this reference 
for further details. 


The Three-Hole Probe 

A drawing of the cobra three -hole probe used in this experiment is shown in Figure A.l. 
The probe axes and nomenclature assigned to the individual openings are also presented. 
The probe was constructed of three silver brazed stainless steel tubes, each having an inside 
diameter of 0.056 cm. The opening of the center tube was normal to its centerline. The 




Figure A.l The cobra three-hole probe 
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openings of the two outer tubes were inclined 30° to their centerline. This arrangement 
resulted in yawing symmetry. 


Yaw-Nulling Calibration 

The main objective of performing a yaw-nulling three-hole probe calibration was to 
empirically determine the relationships between the flow conditions ( p 0 , p) and the three 
measured pressures of the probe (pi, p 2 , Pz)- Since the yaw angle was measured directly, 
the flow total and static pressures were determined as functions of two of the three measured 
probe pressures. To aid calibration, the three pressures measured during probe calibration 
were nondimensionalized as indicated in Equation A.l. Here, p l is the measured pressure of 
port i (i= 1,2,3) and p 0 and p are the local flow total and static pressures. 

C = = fi (Re, M, 7 , a, /?) (A.l) 

Vo ~V 

The pressure coefficients were dependent on Mach number, M, Reynolds number. Re, 
specific heat ratio, 7, the flow pitch angle, a, and the flow yaw angle, /3. The Mach number 
and Reynolds number were interdependent during the calibration and the test, therefore the 
Reynolds number was removed from the list of independent variables. The specific heat 
ratio was also removed since both the calibration and the experimental were performed in 
air. Finally, the experimental flow field was assumed to contain no pitch component and 
the yaw angle was measured directly, so both flow angles were omitted from the list of 
independent variables. This left Mach number as the only calibration variable. 

The calibration was accomplished by varying Mach number while measuring the probe 
pressures. The probe openings were positioned in the center of the jet flow field produced by 
the calibration facility. During calibration, the total pressure was measured by the Pitot probe 
in the calibration facility. Static pressure was simply the test cell atmospheric pressure. The 
probe was nulled to produce equal measured pressures for ports 1 and 2. Typical calibration 
results are shown in Figure A.2. The three coefficients are plotted as functions of total-to- 
static pressure ratio. By varying Mach number, the total-to-static pressure ratio was also 
varied. The coefficient for the center opening, Cpz, was independent of Mach number. The 
coefficients for the side openings , Cp\ and Cpo, were equal (since p\ = p 2 ) and varied 
linearly with total-to-static pressure ratio for the applicable experimental range. 
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A Taylor series expansion was used to approximate the relation between the flow 
conditions and the measured pressures. The approximations for each pressure coefficient 
are given in Equations A .3 and A. 2 . Notice that only constant and linear terms of the 
expansion are included. 

Cp i )2 = b 0 + h (^\ + 0 +... = 6 0 + 6i (^j (A. 2 ) 

Cp3 = Co + + = G ° ( A -3) 

The calibration coefficient ao was simply equal to Cpz, which was one. The calibration 
coefficients 60 and 61 were found using a least squares procedure. This procedure required 
finding the coefficients which minimized the error, given by Equation A. 4 . The minimization 
was accomplished by taking the derivative of the error function (Equation A. 4 ) with respect 
to the coefficients 6 0 and 61 and equating the result to zero. The equations developed from 



Figure A .2 Variation of the three-hole probe 
pressure coefficients with total-to-static pressure ratio 
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this minimization process are presented in Equations A.5 and A.6. Here, N refers to the 
total number of calibration measurements. These two equations were solved simultaneously 
to obtain bo and b\. 


N 

* = £ 


71=1 


Cpi,2 ~ ( ^0 + &1 


P 0 

P 


T 2 


(A.4) 


60JV.t-i.lW7) =E(°Ph 2 ), 

n = 1 V P ' n n = 1 


(A.5) 


N / \ N / \ 2 N / \ 

b ° E ( 7 ) +b 'E ( 7 ) = £«**).(*) 

n= i \ r / n n=1 \ r / n n=1 \ P / 


(A.6) 


Data Reduction 

The local total and static pressures were found from the Taylor series approximations 
to the calibration pressure coefficients. The total pressure was simply the center opening 
pressure, pz, as stated in Equation A.7. The static pressure was found by equating the 
pressure coefficient Cpi ,2 (Equation A.l) to the Taylor series expansion of C'pi ,2 (Equation 
A.2) and then solving explicitly for the static pressure. The resulting quadratic, given in 
Equation A.8 , was formulated. One root of the quadratic was within the calibration region 
while the other root was outside the range of calibration. No difficultly arose in determining 
the correct static pressure. The total and static pressures were then used to determine the 
local Mach number as given by Equation A.9. Since the yaw flow angle varied only slightly 
in the experimental flow, no further data reduction was necessary. 


Po = Pz 


(A.7) 


P = 


— b ± \/& 2 — 4ac 


2a 


: a = 1 - 6 0 , b= (b 0 - bi)pz - pi,2, c = b\p\ (A.8) 




1 




M = 


1/2 


(A.9) 
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APPENDIX B 

FIVE-HOLE PROBE CALIBRATION AND DATA REDUCTION TECHNIQUES 

This appendix contains supplementary information on the calibration methods and data 
reduction techniques for a five-hole probe. Both yaw-nulling and non-yaw-nulling methods 
are reviewed. The information covering three-hole probes can be found in Appendix A. Both 
the three- and five-hole probe calibration techniques used in this thesis closely follow the 
methods developed by Reichert [33]. The reader in referred this reference for further details. 

The Five-Hole Probe 

A drawing of the five-hole probe used in this experiment is shown in Figure B.l. The 
probe axes and nomenclature assigned to the individual openings are also presented. The 
probe was constructed of five silver brazed stainless steel tubes, each having an inside 
diameter of 0.056 cm. The opening of the center tube was normal to its centerline. The 
openings of the four outer tubes were inclined 45° to their centerlines. This arrangement 
resulted in yawing and pitching symmetry. 




Figure B.l Pictorial of the five-hole probe 
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Yaw-Nulling Calibration 

Reichert previously developed an effective yaw-nulling calibration technique for five- 
hole probes [33]. Reichert’s calibration method permits the unknown flow conditions within 
the test flow field to be explicitly calculated during data reduction. The flow conditions 
are determined from measured probe pressures and the calibration coefficients established 
empirically from actual calibration data. This empirical calibration technique is quit different 
than the theoretical calibration approaches used by Vakili et al. [9] in their studies. 
The following section summarizes Reichert’s calibration technique and is included for 
completeness. 

The yaw-nulled five-hole probe calibration empirically determined the relationships 
between the flow conditions (p 0 > p, and a c ) and the five measured pressures of the 
probe (pi, • • • ,ps)- The yaw angle, /?, was measured directly buy nulling the probe in 
the experimental flow field. In practice, the flow conditions (po> p, and a c ) were found from 
three of the five measured pressures. 

The pressures measured during probe calibration were nondimensionalized as indicated 
in Equation B.l. Here, p ( - is the measured pressure of port i (i= 1,..,5) and p 0 and p are 
the local flow total and static pressures. 

Cpi = = fi(M, Re,~f,a c ,0) (B.l) 

Po-P 

As in the three-hole probe calibration, the pressure coefficients were dependent on Mach 
number, M, Reynolds number. Re, specific heat ratio, 7 , measured flow pitch angle, a c , and 
the flow yaw angle, /?. Mach number and Reynolds number were linked during testing and 
calibration, therefore Reynolds number was removed from the list of independent variables. 
Both the calibration and the experiment were performed in air, so the specific heat ratio was 
also removed. Finally, the yaw angle was omitted since the probe was nulled. Therefore, 
Mach number and pitch angle were the independent calibration variables. 

The calibration was accomplished by varying Mach number and pitch angle while 
measuring the probe pressures. The probe openings were positioned in the center of the 
jet flow field produced by the calibration facility. The total pressure was measured by the 
Pitot tube, shown in Figure 113.1. Static pressure was the test-cell atmospheric pressure. The 
probe was nulled to produce equal measured pressures for ports 2 and 4. Typical calibration 
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results are shown in Figure B.2. The five coefficients are plotted as functions of pitch angle, 
a c , for three Mach numbers (M = 0.2, 0.4, 0.6). The pressure coefficients for this probe were 
virtually independent of Mach number in the range of Mach numbers tested. Consequently, 
the pressure coefficients were functions of pitch angle alone. 

Reichert points out that care must be taken when defining the pitch angle. Probe 
pressures responded to the velocity relative to the probe axis, not to the velocity relative 
to the coordinates of the measurement plane, as depicted in Figure B.3. The true pitch 
angle (pitch angle relative to the measurement plane), a, was related to the measured probe 
pitch angle, a c , and the deflection angle, a 0 , as given in Equation B.2. All of these angles 
were measured positive in the clockwise direction. The true pitch angle (a) was defined to 
be the angle between the velocity vector to the measurement plane normal. The measured 
probe pitch angle (a c ) was the angle between the velocity vector to the probe axis and was 
determined from calibration relationships. The deflection angle (a 0 ) was the angle between 
the probe axis to the measurement plane normal. The deflection angle was a result of 



Figure B.2 Variation of pressure coefficients with pitch angle for various Mach numbers 
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imperfect instrumentation manufacturing, occasional bumping of the probe tip during testing 
and/or aerodynamic loading of the probe while cantilevered out into the experimental flow 
field. One simple way of calculating the deflection angle, a 0 , is presented in the Yaw-Nulling 
Data Reduction Section. It is important to remember that the calibration was based on the 
probe pitch angle, a c , and not the true pitch angle, a. 

a = a c + a 0 (B.2) 

Taylor series expansions were used to approximate the relations between the flow 
conditions and the measured pressures. The approximations for each pressure coefficient are 
given in Equations B.3— B.7. Simplifications to the individual calibration coefficients within 
each approximation are included. These simplifications were due to probe symmetry and 
are extensively reviewed in Reference [33]. 


^ 9 

Cpi = a i + a 2 (x c + aza c 

(B.3) 

Cp 2 = ai+ a 4 a 2 c 

(B.4) 

Cpz = aj — a 20 : c + 03 a“ c 

(B.5) 



Figure B.3 Five hole probe pitch angles 
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Cp4 = G] + 040;, 


(B.6) 


Cp5 = 1 + a,5a 2 c (B.7) 

Reichert proved that relationships between these approximations and the flow conditions 
(Pq, p, and a c ) existed. The first, third and fifth expansions were need to develop the 
relationships. Hence, only four of the five calibration coefficients (aj, a2, 03 and a{) were 
determined. These four calibration coefficients were found using a least squares procedure. 
This procedure required finding the values of the coefficients which minimized the error, 
represented by Equation B.8. Here, J is the total number of calibration Mach numbers and 
I\j is the total number of calibration pitch angles for each Mach number. 


E = 


E EE[ C P i(Mj,a Ck ) Cpi(Mj,a Ck ) 

*=1,3,5 j = 1 fc=i 


(B.8) 


The minimization was accomplished by taking the derivative of the error function 
(Equation B.8) with respect to the coefficients (ai, <22, <23 and as) and equating the results 
to zero. The equations developed from this minimization process are presented in Equations 
B.9— B.12. Equations B.10 and B.12 were solved independently to obtain the coefficients 
02 and <25. Equations B.9 and B.ll were solved simultaneously to obtain the coefficients 
ai and 03. 


J J Kj j Kj 

a, Y, Kj + “3 E £ = 5 E £ [Cpi(Mi, a c .) + CMM„ a«)l (B. 9 ) 

j = 1 j = 1 fc=l L j = 1 k= 1 

J Kj J Kj 

“2 E E “I = 2 £ £ - Cn(M v « c JK, (B. 10 ) 

j = 1 *=1 2=1 k=l 


J Kj J Kj J K, 

ai ^2 £ °i + °3 = 2 a c k ) + Cps{Mj, a c J]a4 

j— 1 k=l j= 1 k=l " j= 1 k= 1 


(B.ll) 
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J K, J K, 

+ fl 5EE tt c, 


j = 1 k-1 j = 1 *=1 


7 AO 

^ 1 ^ a Ck ) cx~ k 


j= 1 k = l 


(B.12) 


Yaw-Nulling Data Reduction 

The data reduction procedure determined the local total pressure, static pressure and 
Mach vector components from the Taylor series approximations to the calibration pressure 
coefficients and the measured yaw angle. The reduction process began by determining the 
pitch angle from a ratio of differences of measured pressures, represented by Equation B.13. 
This equation is a quadratic in a c . Although two roots existed for Equation B.13, one root 
was always within while the second root was well outside the range of calibrated pitch 
angles. No difficulty arose in deciding which root was correct. 


— (Pi Ps) _ (B.13) 

2p5-(Pi+P3j (a 5 -a 3 )o;“ + 1 — ai 

The flow local total and static pressures were found using Equation B.14 and the 
previously determined value of a c . These total were then used to determine the magnitude 
of the local Mach number, as given by Equation B.15. 


a\ + a 3 o? c 

1 — ai — a 3 or c 

Vo 


’£(Pi +!*)" 

_ 1 + a 5 

-a s a 2 c 

. P . 


P5 


(B.14) 


M = 



(B.15) 


The last step in the data reduction procedure was to calculate the components of the 
Mach vector relative to the measurement plane. The three components of a unit vector 
representing the flow direction were defined as e z -, e y - and e z - . The component e x - was in 
the direction normal to the measurement plane and in the direction of streamwise flow. The 
components e y - and e z - were projections of the unit vector onto the measurement plane and 
were mutually orthogonal. These three components were found with Equation B.16. 



"1 

0 

0 

= 

0 

COS (f> 

— sin (j> 


0 

sin <f> 

COS (j> 


cos (p + fa) cos (a c + a 0 ) 
sin (p + p 0 ) cos (a c + a 0 ) 
sin (a c + a 0 ) 


(B.16) 
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In Equation B.16, the second matrix determined the flow direction relative to the probe 
coordinate system (Figure B.l). The first matrix transformed the flow components from 
the probe coordinates to the measurement plane coordinates. Two new variables were 
introduced: the reference yaw angle, /3 0 , and the orientation angle, <f>. The reference yaw 
angle, f3 0 , was assigned the negative value of the yaw angle when the probe was nulled at 
the centerline. Therefore, all subsequent yaw angles along a radial traverse were referenced 
to /3 0 . The orientation angle, <f>, was simply the rotation angle about the centerline within 
a measurement plane (Figure III.4). 

Before the three separate components were found, the deflection angle, a 0 , was 
determined. The equation which established the deflection angle is given in Equation B.17. 
This equation was derived by equating the true pitch angle at the centerline when measured 
from two separate orientation angles. Equation B.17 assumes that the true pitch angle at the 
centerline was constant and that the deflection angle did not change while rotating from the 
first to the second orientation angle. Care was taken while handling the probe in order to 
insure the later approximation was a good assumption. 


a 0 = tan 


-l 


sin (<fo 4- Q C2 ) ~ sin (<j>i + q Cl ) 
cos (<f>i + a c , ) - cos (<f > 2 + a C2 ) 


(B.17) 


In practice, the deflection angle was calculated only once as described above. From this 
deflection angle value, the true pitch angle at the centerline, a Ce( , could then be calculated, 
and hence, the left side of Equation B.16 was known at the centerline. For an arbitrary 
orientation angle, therefore. Equation B.l 8 applied at the duct centerline. 


ol 0 = sin 1 (e 2cJ ') - a c - <f> (B.18) 

Equation B.18 was used in general to find the deflection angle along a single radial 
traverse. In Equations B.17 and B.18, the required measure centerline pitch angles were 
previously found from Equation B.13. 


Non- Yaw-Nulling Calibration 

In the remaining sections, the yaw-nulling calibration technique is extended to a non- 
yaw-nulling method. As before, the calibration permits the unknown flow conditions in the 
test flow field to be found from the five measured pressures and the calibration coefficients. 
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However, the non-yaw-nulling data reduction technique requires an iteration scheme to solve 
the non-linear relationships which arise. The following section summarizes details involved 
with a forth-order approximation calibration procedure. 

The non-yaw-nulled five-hole probe calibration empirically determined the relationships 
between the flow conditions (p 0 , p, ct c , and /3 C ) and the five measured pressures of the probe 
(pi, • • • ,ps). The pressures measured during probe calibration were nondimensionalized as 
indicated in Equation B.l. The pressure coefficients were dependent on Mach number, M, 
Reynolds number. Re, specific heat ratio, 7 , the measured flow pitch angle, a c , and the 
measured flow yaw angle, j3 c . As in the yaw-nulled probe calibration, Reynolds number and 
specific heat ratio were removed from the list of independent variables. However, since the 
probe was insensitive to changes in Mach number (determined in the yaw-nulled calibration), 
Mach number was also omitted. Therefore, the measured pitch angle and yaw angle were 
the independent calibration variables. 

The calibration was accomplished by varying pitch angle and yaw angle while measuring 
the probe pressures. The probe openings were positioned in the center of the jet flow field 
produced by the calibration facility. The total pressure was measured by the Pitot tube, shown 
in Figure IELl. Static pressure was the test-cell atmospheric pressure. Typical calibration 
results are shown in Figure B.4. The five coefficients are plotted here as functions of 
measured pitch angle, a c , and measured yaw angle, (3 C , for one Mach number (M = 0.4). 

Similar to the yaw-nulled calibration, care was taken when defining the pitch angle. 
However, defining the yaw angle was less restrictive since the probe was rotated about its 
z-axis for the first measurement at the centerline of the duct (ie... a reference yaw angle, 
f3 0 , was defined as in the yaw-nulled calibration). The true pitch angle, a, was related 
to the measured probe pitch angle, a c , and the deflection angle, a 0 , as given in Equation 
B.2 and shown in Figure B.3. Again, all of these angles were measured positive in the 
clockwise direction. 

Taylor series expansions were used to approximate the relationships between the flow 
conditions and the measured pressures. An approximation for a general pressure coefficient is 
given in Equation B.19. For a forth-order approximation, there exists 15 coefficients for each 
probe hole. Therefore, an “unsymmetric” five-hole probe has 75 coefficients which would 
need to be determined (... a task to which no one would look forward) However, because 


Pressure coefficient, Cp 
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Figure B.4 Variation of pressure coefficients with pitch and yaw angle for 0.4 Mach number 
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the five-hole probe was symmetric, most of the calibration coefficients were eliminated. The 
simplifications to the calibration coefficients, due to probe symmetry, are now reviewed. 


There existed three symmetry conditions which held for the five-hole probe. The first was 
a reflection across the jtz-plane. The second was a reflection across the xy-plane. The third 
was a 90° rotation about the x-axis. These symmetry conditions had important consequences 
for the Taylor series approximations. The simplifications to the calibration coefficients from 
each symmetry condition were developed by Reichert and are summarized below. 

1. Symmetric with respect to reflection across the arz-plane. Therefore 



a i'10 a c + a iO10c+ 

Ci2 0Oc + + a*02/?c + 

fli'30<^c + amalfic + ai\20c c fil 4- g,03/2c + 

a i40 Q c + 0.iz\OL Z c f3 c -|- a t 22a 2 c Pc + a il3 Q c0c + a i04@t + * * * 


(B.19) 


a. Cpi(a c ,0 c ) = Cpi(or c ,-/?c) 


• flioi = 0 

• am = 0 

• «121 = 0 


• G103 = 0 

• ai3i = 0 

• «113 = 0 


b. Cp2(ot c ,/3 c ) — Cp3(o; c , /? c ) 


• 0301 = 0 

• 0311 = 0 

• 0321 = 0 


* 0303 = o 

* O331 = 0 

* 0313 = 0 


c. Cp 5 (a c ,/3 c ) = Cp 5 (a c ,-/3 c ) 


• 0501 = 0 

• 0511 = 0 

• 0521 = 0 


• a 503 = 0 

• 0531 = 0 

• 0513 = 0 


I 
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d. Cp2(oi c ,0c) = Cp 4 (ac,-/?c) 


• 0200 = «400 

’ 0203 = — 0403 

• 0210 = 0410 

0240 = O440 

* «201 — — 0401 

0231 = — 0431 

• 0220 = 0420 

0222 = O422 

* 0211 = —0411 

0213 = -0413 

• 0202 = O402 

• 0221 = — 0421 

0204 = O404 


2. Symmetric with respect to reflection across the xy-plane. Hence 

a. Cp 2 {a c ,Pc) = Cp 2 (-oc c ,Pc) 


• O 2 l0 = 0 

0230 = 0 

* 0211 = 0 

0231 = 0 

• 0212 = 0 

0213 = 0 

b. Cp 4 (a c ,/ 3 c ) = Cp 4 (-a c ,/? c ) 

• 0410 = 0 • 

0430 = 0 

• o 4n = 0 

0431 = 0 

• 0412 = 0 

O413 = 0 

c. Cp 5 ( a c , P c ) = Cp 5 (-a C , 0 C ) 
• 0510 = 0 

' 0530 = 0 

• a 5 n = 0 

' 0531 = 0 

• 0512 = 0 

' O513 = 0 

d. Cpi { a e , Pc ) = Cp 3 (-a c ,/?c) 
• aioo = 0300 

0112 = — 0312 

* Ono = — 0310 

0140 = O340 

• aioi = 0301 

0131 = —0331 

• 0120 = 0 320 

0122 = O322 

• am = —0311 

0113 = —0313 

• aio2 = 0302 

* O130 = —O330 

* O104 = O304 
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3. Symmetric with respect to 90° rotation about the z-axis. Consequently 


a. Cp\{0,f} c ) = Cp 2 (-/3 c ,0) ^ Cpi(a c ,0) = Cp 2 (0,a c ) 


• aioo = <*200 

• ano = 0201 

• <*ioi = “<*210 

• 0120 = <*202 

• 0102 = <*220 


<*130 = <*203 
<*103 = —<*230 
<*140 = <*204 
<*104 = <*240 


b. Cp 5 (0,/3 c ) = Cp5(-/3 c ,0) and Cps(<*c,0) = Cp 5 ( 0 , a c ) 


<*510 = <*501 
<*501 = —<*510 
<*520 = <*502 


<*530 — <*503 
<*503 = —<*530 
<*540 = <*504 


Finally, the aerodynamic stagnation condition implied 


Cp 5 (0, 0) = 1 => 0500 = 1 • 


(B.20) 


The resulting simplifications to the each of the individual coefficients are listed in Table 
B.l. Notice that only 14 independent calibration coefficients were needed for the Taylor 
series approximations. 

Defining the vectors Cp jk , u and a 7 ' (Equations B.21 and B.22), and the matrix My* 
(Equation B.23) enables a system of equations to be written, as represented in Equation B.24. 


CPjk = 


'Cpi (a C; , /?<;*) ' 


'Cpi (a C] ,Pc k )' 


‘O' 

Cp 2 (a Cj ,/3 Ck ) 


Cp 2 (a CJ ,p Ck ) 


0 

Cpz (&cji Pck) 

Cpy* = 

Cp 2 (a C] ,Pc k ) 

u = 

0 

Cpi{oi c ,,fic k ) 


Cp4(a c ,,p ck ) 


0 

_Cp 5 (a Cj ,p Ck ) _ 


-Cp 5 {a Cj ,Pc k ) - 


_ 1 _ 


(B.21) 


<*110 <*120 <*102 


<*502 <*130 


a r =[<*100 
[<*221 


<*140 <*122 <*222 <*104 <*540 


<*112 ] 
<* 522 ] 


(B.22) 
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Table B.l Results of probe symmetry for calibration coefficients 



i = 1 

i = 2 

i = 3 

i = 4 

i - 5 

a i00 

cioo 

ClOO 

ClOO 

ClOO 

1 

0*10 

CllO 

0 

— ano 

0 

0 

0*01 

0 

Clio 

0 

— Clio 

0 

Oi 20 

0120 

Cl02 

Cl20 

C102 

a520 

Gill 

0 

0 

0 

0 

0 

a*02 

C102 

C120 

C102 

C120 

0520 

Oi30 

«130 

0 

-0130 

0 

0 

«i21 

0 

a 221 

0 

-C221 

0 

Ct'12 

C112 

0 

— C112 

0 

0 

Oi03 

0 

ai30 

0 

-ai30 

0 

Ot40 

«140 

ai04 

C140 

C104 

0540 

Ot31 

0 

0 

0 

0 

0 

ai 22 

C122 

C 222 

C122 

a 222 

a 522 

Oil3 

0 

0 

0 

0 

0 

O?04 

aio4 

C140 

Cl 04 

ai40 

0540 


^A-jk — (^Cj ’ ) 


1 

<*c 

<4 

01 

0 

<4 

< *C0c 


1 

0c 

0C 

<4 

0 


01 

0 


1 

-<*C 

<4 

0C 

0 


-<4 

-<x c 0c 


1 

-0c 

0C 

<4 

0 


-01 

0 


0 

0 

0 

0 

<4 + 0c 


0 

0 



0 

4 

O0 c 0l 

0 

0 A c 


0 

0 

eipc 

0t 

0 

a c0c 

<4 


0 

0 


0 

<4 

O0 c 0l 

0 

ft 


0 

0 


C0 c 0 c 

ft 

0 

<4 ft 

<4 


0 

0 


0 

0 

0 

0 

0 

a c 

+ 0c <4 Pi 


(B.23) 


/N 

Cp ]k « Cp jk = u + Mj k a 


(B.24) 
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Relationships between the Taylor series approximations to the pressure coefficients 
and the flow conditions (po> P, ct c , and /3 C ) existed. All five expansions were needed to 
develop these relationships and therefore all 14 calibration coefficients were required to be 
found. These calibration coefficients were again determined using a least squares procedure. 
The procedure required finding the values of the coefficients which minimized the error, 
represented by Equation B.25. 


The minimization was accomplished by taking the derivative of the error function 
(Equation B.25) with respect to each of the calibration coefficients and equating the results 
to zero (Equation B.26). The system of equations developed from this minimization process 
is presented in Equation B.27. Here the vector Cp^- is given in Equation B.21. Equation 
B.27 was solved to obtain the calibration coefficients in vector a. 


The data reduction procedure determined the local total pressure, static pressure and 
Mach vector components from the Taylor series approximations to the calibration pressure 
coefficients. The reduction process began by determining the pitch and yaw angles from a 
ratio of differences of measured pressures, represented by Equations B.28 and B.29. Here 
D is given in Equation B.30 



(B.25) 


»=1 j k 




(B.27) 


Non- Yaw-Nulling Data Reduction 


4ps - (pi + P 2 + Pz + Pa) 


2(pi Pz) 


Quo + Q130 <*c + Q 112 Ocfil 

D 


(B.28) 


4p5 - (pi + P2 + P3 + Pa) 


2(P2 - Pa) 


Qlioffc + Ql30/?c + Q221 ttcffc 

D 


(B.29) 
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D = 1 — <*100 + ^<*520 — ^(<*120 + <* 102 H ( Q c + Pc) 

+ ^<*540 — ^(<*140 + <* 104 )^ ( Q c + Pc) (B.30) 

+ ^<*522 — ^( a 122 + 0,222)) (&cPc) 

The right-hand-sides of Equations B.28 and B.29 are fourth-order approximations to the 
respective / and g functions. These equations are non-linear and coupled. An initial guess 
to the measure pitch and yaw angles was obtained by developing and solving second-order 
approximations to the / and g functions. The second-order approximations are given in 
Equations B.31 and B.32. 


2(pi -P3) 

4p5 - (Pi + P2 + P3 + P4) 


<*110 Qc 

1 - O 100 + (<*520 - ^(<*120 + <*102)) (<*c + Pc) 


(B.31) 


2(P2 — P 4 ) 

4p5 - (pi + P2 + P3 + P4) 


OUoPc 

1 - <*100 + (<*520 — |(<*120 + <*102))(<*c + Pc) 


(B.32) 


The right-hand-sides of the second-order approximations to the respective / and g 
functions are also non-linear and coupled; however, the non-linear and coupling term 
[op. + Pc) * n the denominators of Equations B.31 and B.32 can be found by solving the 
quadratic given in Equation B.33. Consequently, a first guess to measure pitch and yaw 
angles can be solved for explicitly from Equations B.31 and B.32. 


<*( Q c + Pc) + K a c + Pc) + C — 0 


<* = (/ 2 + g 2 ) ^<*520 — 2 [ a 120 + <* 102 ]^ 

b = 2 (f 2 + g 2 ) ^<*520 — ^[<*120 + <* 102 ]^ (1 + <* 100 ) - <*no 
c= (f 2 + g 2 )( 1 + <*ioo)“ 


(B.33) 
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Non-linear terms in Equations B.28 and B.29 were linearized. The initial approximations 
to the measure pitch and yaw angles were used as a starting point for solving the fourth- 
order approximations represented by B.28 and B.29. The measure pitch and yaw angles 
were obtained iteratively and were determined upon convergence. 

The local total and static pressures were found using Equation B.35 and the previously 
determined value of a c and f3 c . Here, the pressure coefficients, Cp, are defined by Equation 
B.24. The total and static pressures were then used to determine the magnitude of the local 
Mach number, as given by Equation B.15. 


The last step in the non-yaw-nulling data reduction procedure was to calculate the 
components of the Mach vector relative to the measurement plane. This step was completed 
in the same fashion as previously described in Yaw-Nulling Data Reduction and is therefore 
not presented here. 


Cpi234 1 — Cpi234 f Po 


Cp 5 1 - Cps J [ p 


\{Pl +P2+P3 +P4 ) 
P 5 


(B.35) 
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